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We apply the approximate dynamics derived from the Gaussian time-dependent variational principle to the Hamiltonian 
H = \{p x + py) + r;X 2 y 2 , which is strongly chaotic in the classical limit. We are able to calculate, essentially analytically, 
low-lying eigenvalues for this system. These approximate eigenvalues agree within a few percent with the numerical results. We 
believe that this is the first example of the use of TDVP dynamics to compute individual eigenvalues in a non-trivial system 
and one of the few such computations in a chaotic system by any method. 



I. INTRODUCTION 

Gaussian approximations to quantum mechanics have 
been utilized successfully in many contexts including 
quantum field theory jlj], the dynamics of hydrogen 
plasma §], semiclassical propagation methods HQ, 
quantum control [^) and the study of "quantum chaos" 
|^| . The primary motivation for their popularity is sim- 
plicity of computation: Gaussians are easily parame- 
terized by the c number variables specifying the cen- 
troid (average variables) and spread (fluctuation vari- 
ables) and their dynamics are essentially classical, apart 
from the computation of a phase which is a crucial ele- 
ment. Further, Gaussians arise naturally in the coherent 
state representation of quantum mechanics and in the 
N = oo limit (where N is number of degrees of freedom) 
of many-body systems || . There are a variety of Gaus- 
sian approximations including: (1) a variational approxi- 
mation usually derived through the time dependent vari- 
ational principle (TDVP) (2) a recently intro- 
duced quadratic-order Gaussian approximation Q, (3) 
Heller's method which is a non-self-consistent trun- 
cation of (2), and (4) the multiple classical trajectory 
version of Heller's method pl| . A further level of approx- 
imation yields the Gaussian Effective Potential method, 
which consists of an adiabatic elimination of the time- 
dependence of the fluctuation variables in the TDVP dy- 
namics |p^| . 

The wide-spread use of these methods raises the ques- 
tion of their validity and range of applicability |p|. It 
has been argued that in the presence of chaos, semiclassi- 
cal approximations to quantum mechanics should break 
down on a logarithmic time-scale t c ~ ylog(^), where 
A is the largest Lyapunov exponent of the underlying 
classical mechanics, and h is Planck's constant. Com- 
putations with the multiple trajectory Gaussian approx- 



imation has demonstrated that this may be a pessimistic 
estimate @. Recently, the validity of the TDVP approx- 
imation has also been considered. Apart from the chaos 
in the underlying classical system, it has been shown 
that the approximate quantum dynamics derived from 
the TDVP may be chaotic even when the classical limit 
is not This has led to the argument |ll| that the 

TDVP Gaussian approximation fails in the presence of 
chaos. However, it has been shown that this anoma- 
lous chaotic behavior persists even when exact numerical 
computations are made p^| . It is argued, in fact, that 
this chaos is a signature of the complicated nature of 
the spectrum involved in the exact quantum dynamics 
0,0. Further, recent work by Habib Q shows that all 
Gaussian approximations to Schrodinger's equation are 
identical to the same approximation to the classical Li- 
ouville equation, although the classical versions do not 
have any phase information. This result clarifies that fi 
is a kinematical constant in these approximations, pro- 
viding a scale for the "smoothing" of the dynamics and 
reinterprets the "quantum effects" included in Gaussian 
approximations. It does not invalidate the results of the 
quantum Gaussian approximations, although it does re- 
quire these results to be understood in kinematical rather 
than dynamical terms. The result also emphasizes that 
quantum dynamics are better approximated by classical 
Liouville dynamics rather than Hamilton's equations for 
point trajectories JnJ. 

It is thus clear that Gaussian approximations should be 
used and interpreted with caution. However, as we shall 
demonstrate in this paper, the TDVP Gaussian approx- 
imation does yield accurate results even in the presence 
of chaos, in a system where other approximations fail. 
The TDVP dynamics can be used to compute eigenval- 
ues [pOP] through an extension of the Einstcin-Brillouin- 
Keller quantization method pi]] . We use this method 
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to compute eigenvalues for the two-dimensional coupled 
quartic oscillator 
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This system is highly chaotic classically p2|; it was be- 
lieved till recently to be ergodic, and the integrable re- 
gions of phase space occupy less than 0.005% of the vol- 
ume. The quantal Hamiltonian also resists numerical 
analysis; large basis-sets do not suffice for quantization 
p3[ . It is usual for numerical ease to add a term like 
(3(x A + y 4 ) to the potential and to analyze the system in 
the limit j3 —> 0. The traditional methods for semiclas- 
sical quantization fail for this particular system. Firstly, 
the EBK method cannot be applied because of the chaos 
in the system - there are no stable torii to use for quan- 
tization. The alternative methods developed in the field 
of quantum chaos || also do not work. The most promi- 
nent of these, the Trace Formula method, starts from 
Feynman's path-integral representation of quantum me- 
chanics, and through a sequence of stationary phase ap- 
proximations, derives the eigenvalues of the quantum sys- 
tem as poles in a weighted sum over the unstable peri- 
odic orbits of the classical system p3. Recent successes 
in other systems |^4|,^5| notwithstanding, the systematic 
enumeration of the unstable periodic orbits of the classi- 
cal system has not been achieved for this particular sys- 
tem. Another alternative, Heller's approach ||, which 
computes an approximate time-dependent wave function 
{9(t)\ and Fourier transforms the overlap (*(t)|*(0)), 
fails because of the inherent instability of the trunca- 
tion, as argued in a previous qualitative analysis in 
model potentials. In the same paper, we demonstrated 
that the TDVP dynamics restores stability to systems 
where the truncated Gaussian dynamics fail; in fact, the 
TDVP dynamics can be stable even when the classical 
dynamics are unstable. This is precisely what happens 
in the system given by Eq. (Q): An infinite set of un- 
stable periodic orbits is stabilized by the "quantal fluc- 
tuation" terms in the TDVP method. These orbits can 
then be used to compute approximate eigenvalues for a 
symmetry subspacc of this system. Remarkably, these 
are obtained analytically (barring a numerical integral) . 
Further, these results are obtained with = 0, i.e. the 
regime where even large basis-set calculations fail to con- 
verge. The high degree of chaos in the system and the 
fact that we use h = 1 in our calculations suggests the 
naive perspective that the classical-like TDVP approach 
is then far from its region of validity. However, our results 
are extremely accurate when compared with the "exact" 
numerical results over a substantial range for the lowest- 
lying eigenvalues. Thus, while there is no suggestion that 
this method can be always used in the presence of chaos 
to successfully approximate quantum dynamics, our re- 
sults indicate that it can certainly be used with care in 
some circumstances. 



In the next section, we shall briefly review the dynam- 
ical equations for the TDVP method, including the con- 
struction of a quantization rule. In the third section, we 
apply this method to Eq. ([!]), comparing the method in 
the process to the usual semiclassical methods. We then 
discuss the results and argue that the regime of validity 
of the method is the low quantum-number regime, con- 
trary to the usual intuition derived from the Correspon- 
dence Principle that classical-like approximations work 
best in the high quantum number regime ||. We thus 
suggest that the TDVP Gaussian approximation works 
best as a technique complementary to the usual semi- 
classical methods. 



II. TDVP GAUSSIAN DYNAMICS 

The usual derivation of these dynamics proceeds from 
Dirac's Time-Dependent Variational Principle (TDVP) 
^[Tof; this posits an action of the form T = 
j dt(^,t\ih-gi — H\^f,t). The general requirement that 
5T = yields the Schrodinger equation and its complex 
conjugate. The true solution is approximated by restrict- 
ing |\&(t)) to a subspace of the full Hilbert space and 
setting ST — within this subspace. In particular, this 
restriction may be to the space of general Gaussians pj. 
We have derived these same dynamics Jl4|,(§] in a some- 
what more intuitive fashion from Ehrcnfest's Theorem. 
We start with the equations for the centroid variables 
and make a Taylor expansion around the centroid with 
the higher moments of the wave function. These mo- 
ments follow the usual Heisenberg dynamics, yielding in 
general an infinite system of equations. We render this 
system finite by projecting onto the space of Gaussians; 
this system of equations are identical to those derived 
from the TDVP. 

We have been able to represent these dynamics as an 
extended classical gradient system for the average and 
fluctuation variables with dynamical equations 
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where the subscript ext indicates the "extended" poten- 
tial and Hamiltonian and = d n V/dx n \ <x> . The co- 
ordinate variables for the extended Hamiltonians are x, p 
and their canonically conjugate momenta are p, II and are 
related to the moments of the Gaussian ^{x^p, p, II, t) as 
follows: 



(x) = x, (8) 

m^p, (9) 

(AxAp + ApAx) = 2pH, (10) 

{Ax 2m+1 ) = 0, (12) 

pWH^ + p 2 n 2 . (is) 



The first 3 relationships Eqs. (|| - |T^) are definitions and 
Eqs. ([Tl| - |l2| ) are a consequence of the Gaussian ansatz. 
Eq. (|13|) is a kinematical constraint arising entirely from 
choosing the equality in the Uncertainty Principle rela- 
tionship; it is the only way in which h enters this approx- 
imation. 

In this method, the fluctuation and average variables 
are treated on the same footing and the phase space is 
dimensionally consistent: p has the dimensions of length 
and II that of momentum. The geometry of the space 
is thus identical to that of an ordinary classical system 
- it is a Cartesian space, a manifold 1Z(2N) defining the 
extended phase-space. There is also an equation for the 
phase of the wave function: If we define A as |<&, i) = 

exp (j^^-j 1^, t) it is simple to derive from Schrodinger's 
equation the equations for A = Xd + ^g with the first part 

\d = - I dr (H) = -t H ext (14) 
Jo 

corresponding to the dynamical phase. The second part 
is the geometrical phase: 

A G = J dr (ih-^-) = J dr ( PU ~ UP + px). (15) 

For cyclic evolution this is the Aharanov-Anandan form 
of "Berry's phase" ^(|; it depends only on the geome- 
try of the evolutionary path in phase space and can be 
written as 

A G (C)= £ P-dQ, (16) 

where P = (p, II) and Q = (x,p). The equation for 
the phase, along with the Hamiltonian equations of mo- 
tion for the evolution of the wave function parameters 
constitute the TDVP dynamics. This lies on a space 
1Z(2N) x S(l); for the case just considered, TV = 2, the 
result is completely general, however. 



We now provide a constructive argument |t| for obtain- 
ing eigenfunctions and eigenvalues, which is equivalent 
to imposing a single-valuedness constraint on stationary 
wavefunctions (2(]] in the extended phase-space. Note 
that an eigenfunction for the extended dynamics is one 
whose parameters are invariant under the evolution. We 
see readily that a periodic orbit (PO) solution to Hamil- 
ton's equations Eqs. (§-||) is invariant on the 1Z(2N) sub- 
space; however, each point along the PO acquires a phase 
factor during the evolution. The dynamical phase is the 
same for all the points along the PO and can be factored 
as a global phase. The geometrical part Ag for the PO 
is crucial: We note that a PO x Ag such that the peri- 
odic evolution of Ag on 5(1) is commensurate with that 
of the PO on 1Z(2N) is a function invariant on entire 
space 1Z(2N) x 5(1) and is hence an eigenfunction. The 
commensurability of the phase translates to the relation 

Ag(PO) = ^-/ PdQ = n?i, (17) 
2tt Jpo 

where we have used Eq. (|l6|). Thus, the eigenfunction 
is a weighted sum (the weight factor at each point be- 
ing the appropriate geometrical phase) over the points 
of the commensurate periodic orbit and the eigenvalue 
is H ext for that PO. This rule is the same as the "old" 
quantization rule of Bohr and Sommerfeld; however, it 
applies in the extended phase space, as opposed to the 
classical phase space and thus does not have the same 
meaning. In particular, there are no Maslov-Morse cor- 
rections |2l[] to this rule, since there are no singularities 
in the Gaussian representation. It can be shown [p7p|l 
that the "spread" variables p, n explicitly take care of 
these corrections. In general, this quantization condition 
will give results different from the EBK rule (the POs 
are in the extended space) but always incorporates the 
Maslov correction. The extension of this argument from 
POs to invariant torii goes through easily [£0| and leads 
to a general quantization rule: 

-!- I P dQ = mh (18) 

where the closed integral is now taken over the ith irre- 
ducible contour around the torus and the quantum num- 
bers rii are labeled accordingly. This is exactly Einstein's 
generalization [p8| of the Bohr-Sommerfeld rule to invari- 
ant torii. 

The system of equations derived by Heller || for the 
semiclassical evolution of Gaussian wavepackets obtain as 
truncations of equations (|^) and (||) to O(p ) and of equa- 
tions ([|) and (^) to 0(p 1 ). This arguably Q] inconsistent 
semiclassical system of equations destroys the Hamilto- 
nian structure of the dynamics, leading to non-unitary 
evolution (^9). A consistent truncation to 0(p 1 )) for this 
system retains the Hamiltonian structure of the TDVP 
method and has been termed extended semiclassical dy- 
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namics Q : All the advantages of the TDVP method ap- 
plies to the extended semiclassical method, including the 
definition of a Poisson bracket, and the existence of a uni- 
tary propagator and an analytic quantization method. 
Unlike the TDVP method, Heller's Gaussian dynamics 
and the extended semiclassical dynamics arise as "con- 
trolled" first-order expansions; their validity can thus be 
formally evaluated jl3| and these are hence attractive 
approximations. However, the truncation of the dynam- 
ical equations at the term involving the third derivative 
of the potential induces an unphysical instability, where 
the fluctuation variables grow without bound |Ij even 
for simple one-dimensional anharmonic potentials like 
V(x) — x A . In the TDVP dynamics, all orders of deriva- 
tives are maintained with a re-summation of the moment 
expansion under a Gaussian ansatz; this yields behavior 
that is qualitatively similar to the exact long-term quan- 
tal behavior, in particular reproducing the appropriate 
stability. Thus, the TDVP method can stabilize unsta- 
ble periodic orbits, which may then be used as above to 
obtain eigenvalues and eigenfunctions, as we now demon- 
strate. 

III. EIGENVALUES FOR A CHAOTIC 
HAMILTONIAN 

We now turn to the computation of eigenvalues for the 
Hamiltonian 

H = pkt% + l s> 2f +m * + f ) , (19) 

Extensive numerical work (2^]2^] shows that the classical 
limit (O — > O for all operators) is a very strongly chaotic 
system, with few stable periodic orbits in the limit f3 — ► 0. 
It is easy to verify, however, the existence of an infinity 
of unstable periodic orbits along the diagonals of the po- 
tential. We also note here that the Hamiltonian displays 
a simple scaling relationship in energy |23| : A trajectory 
(xi(t), pi(t)) at an energy E\ is related to a trajectory 
(x 2 (f),p2(t)) at an energy E 2 by 

Mr) = (f¥*i(u (20) 

Pa(r) = (|Mpi(*) (21) 

where r is the rescaled time r = (-fj-) - '- This means 
that there is the same degree of strong chaos at all finite 
energies: There is no "transition to chaos". Eckhardt, 
Hose and Pollack have done a careful numerical analysis 
of the quantum system to show the presence of "scars" in 
the eigenfunctions p3| . They state that harmonic oscil- 
lator basis-set quantization with matrices of dimension 
3240 do not provide converging eigenvalues for (3 = 0; 
they have hence used (3 = 0.01 for their analysis. The 



eigenfunctions of this Hamiltonian belong to the symme- 
try classes of the C± v symmetry group which has eight 
elements (four reflections in the axes and diagonals and 
four rotations by 7r/2). The irreducible representations 
of this group split into four one-dimensional represen- 
tations and one two-dimensional representation. They 
have restricted themselves to the four one-dimensional 
representations corresponding to wave functions which 
are A) symmetric under x — > y, x — > —x, B) antisym- 
metric, symmetric, C) symmetric, antisymmetric and D) 
antisymmetric, antisymmetric respectively, and have nu- 
merically obtained low-lying eigenvalues and eigenstates 
for this system. 

We have applied the Gaussian wavepacket methods de- 
tailed above to this system with /3 = 0. Of the three 
methods, the truncated Gaussian methods (both Heller's 
dynamics and the extended semiclassical system) fail, 
yielding unstable motion where the wavepackets spread 
without bound, irrespective of the value of (3. This is eas- 
ily established by noting that there exist one-dimensional 
projections in which this two-dimensional potential re- 
duces to the anharmonic quartic potential considered 
above and the existence of a single unstable direction for 
the dynamics corresponds to instability in the global mo- 
tion. This exposes one particular frailty of the truncated 
Gaussian approximations: they work well in systems that 
are close to harmonic only in the particular sense of being 
potentials of the form V(x) — x 2 +f(x) with the function 
f(x) containing higher polynomials. 

On the other hand, the TDVP method works excel- 
lently for this system; firstly, the dynamics are completely 
bounded. It is an interesting feature of this approach 
that even though there exist classically unbounded or- 
bits along x — or y — (which are precisely what make 
the numerical quantal analysis through basis sets so dif- 
ficult) the inclusion of the re-summed moment terms via 
the variational approach restores stability to the prob- 
lem. This effect has been termed "quantum resuscita- 
tion" in the context of the Gaussian Effective Potential 
The extended Hamiltonian for the TDVP (using 
* = $(a;)$(y)) is 

H ext = ±(pl+ P 2 y + Til + nl) 

+ ^2+^ + \(* 2 +pW + pD- (22) 
8p 2 x %pl 2 

Note that the use of the factored wave functions explicitly 
restricts us to the one-dimensional representations. The 
dynamics of this extended Hamiltonian are, in general, 
chaotic. However, if we consider the subspaces of the one- 
dimensional representation noted above, we find that the 
first and third subspace can be studied by the symmetry- 
reduced version of H ext : 

H=\(p 2 +^) + ^- 2 +\{z 2 +p 2 f (23) 
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where (z,p) and (ft, II) are the canonically conjugate 
pairs. 

We now demonstrate that this symmetry-reduced ver- 
sion of the extended Hamiltonian is explicitly integrable. 
To do so, we make the change of variables to spherical 
coordinates R, 8 defined in the plane by: z — R cos 9; p = 
R sin 9 . This transforms the Hamiltonian to 
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Pi 
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8 sin^ 9 



(24) 



Since the factor multiplying is solely a function of 9, 
this is now in the right form to use Hamilton- Jacobi the- 
ory (|(]]. Through the separability just noted, therefore, 
we introduce Hamilton's characteristic functions Wr and 
Wg and get the Hamiltonian- Jacobi equations: 



<W 



1 



! sin 2 9 



= k 



dWn 
dR 



4 R 2 



(25) 
(26) 



where E (the energy) and k are the constants of separa- 
tion. We form action variables as usual: 
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The 9 integral yields {u| 



Je 



a- 1 



(27) 

(28) 
(29) 
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where a is introduced for convenience through k = 
The R integral is a complicated elliptic integral that can 
be evaluated; however, it cannot be analytically inverted 
to yield the quantization condition. We hence leave it in 
quadrature: 
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dR\ 2E 
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R 4 

T" (32) 



The existence of this integral demonstrates the integra- 
bility of the chosen symmetry subspace of the extended 
Hamiltonian. Since we now have a set of stable invari- 
ant torii in the extended space, we can proceed with 
the quantization as detailed above in a straightforward 
fashion. To wit: Eigenvalues correspond to torii with 
quantized actions in both variables R and 9. We follow 
this prescription by setting Jg equal to a half-integer in 
the above equation [the symmetry-reduced form of the 



Hamiltonian only accumulates half the phase of the ac- 
tual Hamiltonian which is why we use half-integers] ; this 
yields a one-degree of freedom dynamical system in R 
which has only closed orbits. We then proceed as follows: 
We take various initial conditions and numerically inte- 
grate their dynamics over the closed orbit to compute the 
action <f P^dR. The orbits for which the action equals a 
half-integer then correspond to eigenfunctions and their 
conserved energy the associated eigenvalue. We show 
the results for the first 67 eigenvalues in Fig. |J com- 
pared with the numerical results (3^] of Eckhardt et al. 
We now note that the eigenvalues we have calculated es- 
sentially analytically agree within a few percent with the 
numerical results P3| over the substantial range of our 
calculations. 

We emphasize that there is no possibility that the va- 
lidity of the results can be attributed to the minute re- 
gions of stability of the classical phase-space |Q . A mo- 
ment's consideration shows that the use of the factored 
wave functions and the restriction y(t) = x(t) = z(t) for 
the symmetric dynamics corresponds to restricting our 
attention to the diagonals of the potential. That is, the 
dynamics of the Gaussian are restricted such that the 
centroid always travels always along the unstable peri- 
odic orbits along these diagonals. The few classical stable 
periodic orbits that do exist in this system lie far from 
this region and our results cannot be understood as af- 
fected by the presence of these orbits; the wavepacket is 
not influenced by them. 

Further, we note that there is a superficial similarity of 
these results to other work Q] that demonstrates the effi- 
cacy of Gaussian approximations in computing low-lying 
spectral features. However, those results depend on per- 
turbations of a classically integrable Hamiltonians. As 
such, they were able to use standard quantal perturba- 
tion theory. This is not possible for Eq. ([j]). Further, 
since the classical dynamics is ab initio strongly chaotic, 
even the Gustavson-Birkoff quantization technique [p5[ , 
which is an adaptation of classical perturbation theory, 
cannot be applied. In either case, there is nothing to 
perturb around for this system. This also emphasizes 
the non-perturbative aspect of the variational approxi- 
mation. 

It is clear that the TDVP method in this system takes 
advantage of the interesting feature of being able to use 
the infinity of periodic orbits along the diagonals of the 
potential. These orbits are classically unstable, and the 
formal application of the WKB quantization method to 
this unstable periodic orbit yields metastable states |36| , 
where the eigenvalues have a real part (corresponding 
to the energy) and an imaginary part (corresponding to 
the life-time of the state). Apart from the unphysical 
metastability of these states, the approximate energies 
thus obtained are valid only for the first few states - 
we show the limited accuracy in Fig. 1, where we have 
plotted the real parts of the first 10 eigenvalues from 
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this method. This same instability of the periodic orbit 
causes the breakdown of the truncated Gaussian semi- 
classical methods as well. However, the TDVP Gaussian 
ansatz makes dynamics along these periodic orbits sta- 
ble - an example of "quantum resuscitation" - and our 
generalized quantization method succeeds. 

There are more general considerations also: Note that, 
despite the high degree of chaos in the system, the poten- 
tial is relatively benign for Gaussian approximations: It 
has only one minimum and no maxima. However, since 
it is a quartic well, this does not benefit the truncated 
approximations and we must turn to the fully re-summed 
TDVP approximation to take advantage of this structure. 
Second, the high value of h is an advantage in this con- 
text. As recent work |37j demonstrates, at higher values 
of U the details of the classical phase-space structure are 
smoothed out in the quantal dynamics (see figure (2a) of 
|57j), as reflected by the quantal eigenfunctions. Since 
Ti also sets the kinematical scale of smoothness for the 
TDVP Gaussian method |is|] , this means that both the 
exact quantal dynamics and our approximate version are 
effectively occurring in a smoother potential well than 
the classical point dynamics. Both these factors imply 
that distorted Gaussians can be expected to evolve with- 
out excessive error under such circumstances. Further, 
the smoothness permits the weighted superposition of 
Gaussians inherent in the TDVP eigenfunction ansatz to 
yield accurate results: The detailed structure of the wave 
functions may be argued to contribute rapidly oscillating 
terms that affect the computation of averages of observ- 
ables in these states in a small way. Thus, eigenvalues, 
which are the average of the Hamiltonian operator, may 
be quite accurate even when the detailed dynamics of the 
variational approximation are not so. This also clarifies 
that the method is expected to be valid in the regime of 
the lowest-lying eigenvalues, contrary to the usual Corre- 
spondence Principle regime where classical-like approxi- 
mations work best in the high quantum number regime 
. This is supported by our results: they deviate slowly 
away from the numerical results as the quantum number 
increases. 

Our results thus show that the degree of chaos in a clas- 
sical problem does not necessarily limit the validity of the 
TDVP Gaussian dynamics. In a general case, chaos in 
the classical dynamics of a system coexists with compli- 
cated potential energy surfaces; Gaussian methods would 
hence apply only over a small range of parameters. How- 
ever, there are systems that may be extremely chaotic, 
but possess the appropriate structure that enables the 
TDVP Gaussian approximation to work over a much 
larger range, especially in the high Ti regime; this explains 
the many successful applications for instance. It is 
clear that these arguments need to be explored carefully 
in more situations. 

In summary, we have analytically (barring a numeri- 
cal integral) computed the lowest-lying eigenvalues of a 



classically strongly chaotic system. We believe that this 
is the first example of the use of TDVP dynamics to 
compute individual eigenvalues in a non-trivial system 
and one of the few such computations in a chaotic sys- 
tem by any method. These results compare extremely 
favorably with numerical results and show that the lim- 
its of validity of the TDVP method are not necessarily 
set by the degree of chaos in the classical system. Ac- 
curate quantum dynamical simulations are quite difficult 
and there is a great need for valid approximations p8| . 
Heller's truncated Gaussian approximation has already 
been shown || to be extremely useful for systems that 
can be explicitly written as perturbations around a har- 
monic minimum. It is intuitive that Gaussian approxi- 
mations should continue have validity in potential wells, 
even in the presence of chaos; however, the truncated 
Gaussian methods fail in anharmonic systems. We be- 
lieve that the TDVP Gaussian method is an excellent 
candidate for approximate quantum dynamics in these 
and other situations. 
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FIG. 1. Comparison of lowest eigenvalues 
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